Magnetic moment non-conservation in magnetohydrodynamic turbulence models 
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The fundamental assumptions of the adiabatic theory do not apply in presence of sharp field 
gradients as well as in presence of well developed magnetohydrodynamic turbulence. For this reason 
in such conditions the magnetic moment ^ is no longer expected to be constant. This can influence 
particle acceleration and have considerable implications in many astrophysical problems. 

Starting with the resonant interaction between ions and a single parallel propagating electro- 
magnetic wave, we derive expressions for the magnetic moment trapping width (defined as the 
half peak-to-peak difference in the particle magnetic moment) and the bounce frequency loi,. We 
perform test-particle simulations to investigate magnetic moment behavior when resonances over- 
lapping occurs and during the interaction of a ring-beam particle distribution with a broad-band 
slab spectrum. 

We find that magnetic moment dynamics is strictly related to pitch angle a for a low level of 



magnetic fiuctuation, SB/Bq = (10 



10" 



where Bo is the constant and uniform background 



magnetic field. Stochasticity arises for intermediate fluctuation values and its effect on pitch angle 
is the isotropization of the distribution function /(a). This is a transient regime during which 
magnetic moment distribution /(/x) exhibits a characteristic one-sided long tail and starts to be 
influenced by the onset of spatial parallel diffusion, i.e., the variance {(Az)^) grows linearly in time 
as in normal difi^usion. With strong fluctuations f{a) isotropizes completely, spatial diffusion sets 
in and f{li) behavior is closely related to the sampling of the varying magnetic field associated with 
that spatial diffusion. 



I. INTRODUCTION 

In this paper we study magnetic moment fj, conserva- 
tion for charged particles in presence of a single electro- 
magnetic wave as well as in presence of turbulent mag- 
netic fields having one dimensional spectra comparable 
to those measured in the solar wind. Magnetic moment 
conservation is an important topic in plasma physics. In- 
deed, some of the most commonly used theories that de- 
scribe particle motion in perturbed magnetic fields are 
based on the assumption that particles magnetic moment 
is on average constant over a gyroperiod. When /i is not 
conserved, this approximation is not allowed and its ef- 
fects can have a bearing on several astrophysical phenom- 
ena such as coronal heating, cosmic ray transport, tem- 
perature anisotropics observed in the solar wind [l| and 
particle acceleration near reconnection sites . Further- 
more this issue is strictly related to particle confinement 
in plasma machines and dynamically chaotic systems Q . 
Therefore we want to established the validity range of the 
adiabatic approximation and the key mechanisms that 
regulate magnetic moment non-conservation. 

The guiding center approximation (3] splits particle 
motion into the motion of the guiding center and the 
gyromotion around it. When analyzing charged particle 
motion in nonuniform electromagnetic fields, we would 
like to neglect the rapid and relatively uninteresting gy- 
romotion, focusing instead on the far slower motion of 
the guiding center. Averaging the particle equation of 
motion over the gyrophase, we obtain a reduced equa- 
tion that describes the guiding center motion. In the 
non-relativistic case the equation of motion of the guid- 



ing center in the direction parallel to the magnetic field 
reads 



dt 



-IT = -A<V||B + (7-B| 



'II' 



(1) 



where particle magnetic moment is defined as fi = v\/ B 
and V|| = (B • V) is the spatial derivative along the 
field direction. In the perpendicular direction the guiding 
center drifts with the velocity 



VD 



F X B 



(2) 



where F = [gE — ^VB — (m?;|)V||B] is the total force 
acting on the guiding center, averaged over a gyroperiod, 
in the (non-inertial) frame co-moving with the guiding 
center. Therefore, as long as a particle moves through 
slowly varying electric and magnetic fields, its guiding 
center behaves like a particle with a magnetic moment /i 
conserved. 

This approximation is valid when the smallest length- 
scales of the electromagnetic fields are much larger than 
the particles Larmor radius, i.e., when particle magnetic 
moment is a constant of motion on average over the par- 
ticle gyroperiod. This corresponds to the well-known 
Born-Oppenheimer approximation in quantum mechan- 
ics. This description for particle motion in a non-uniform 
magnetic field is also useful for numerical simulations. 
Indeed direct simulations of kinetic equations (Vlasov, 
Boltzmann) with a large magnetic field require the nu- 
merical resolution of small spatial and time scales in- 
duced by the gyration along the magnetic field. The 
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guiding center approximation, as well as gyrokinetics, are 
approximate models describing particle motion in pres- 
ence of a strong magnetic field. However, the assumption 
that the scale of variation of the magnetic field is much 
larger than the particle Larmor radius can break down 
in presence of turbulence. Turbulent magnetic fluctua- 
tions are observed in space plasmas in practically all en- 
vironments and at all scales. Furthermore the presence 
of waves in collisionless plasma introduces through wave- 
particle interactions a finite dissipation. In this case it 
seems invalid to resort to a guiding center theory. 

When the amplitude of the magnetic fiuctuations is 
lower than that of the mean magnetic field (averaged over 
the fluctuations time-scale), a perturbation approach 
called quasilinear approximation is applicable ^tJ\. In 
this case the resonant fluctuations make the dominant 
contribution to particle scattering. The resonance condi- 
tion for wave-particle interaction is given by: 

uj — fc||W|| = nVL (3) 

where lj is the wave frequency, fc|| and U|| are respectively 
the wavevector and the particle velocity along the mean 
magnetic held Bo, and — qB/ni is the particle gyrofre- 
quency. Landau resonance [3] is found at n = 0, while 
n = ±1, ±2, ... are the cyclotron resonances. In linear 
theory these resonances are represented by delta func- 
tions. In presence of well-developed magnetohydrody- 
namic turbulence we expect that the discrete resonances 
to be significantly broadened due to the rapid decorrela- 
tion of the waves phases in strong turbulence . 

The particle reaction to the perturbation is always pe- 
riodic except when condition ^ is satisfied. In this case 
the perpendicular electric force due to the wave remains 
in phase with the particle cyclotron motion and particle 
reaction is secular or resonant and, over short times, non- 
oscillatory. The secular electric force acting on a given 
particle is constant over a particle gyroperiod, so that 
the magnetic moment is no longer conserved. 

Charged particles are scattered by their interaction 
with the waves and undergo pitch angle diffusion. The 
pitch angle, 9 = arctan(?;j_/w||), is the angle between the 
direction of the magnetic field and the particle's helical 
trajectory. Scattering from magnetic fluctuations causes 
the distribution of pitch angle cosine, a = v\\l\v\, to be- 
come isotropic. Magnetic moment, /i, is formally related 
to the time averages of the cosine of pitch angle by: 



We therefore expect the behavior of magnetic moment to 
be strongly related to pitch angle behavior. 

II. STOCHASTIC MOTION, TRAPPING 
WIDTH AND RESONANCE OVERLAPPING 

Wave-particle interactions usually involve multiple res- 
onances. Particle motion is substantially different de- 



pending on when these resonances overlap or not. Nu- 
merical simulations show a complex behavior that can- 
not be approached analytically, e.g., it is not possible to 
write an equation for the evolution of particles distribu- 
tions when two resonances overlap ^lOj . Such motions 
in the presence of overlapping resonances are commonly 
labeled stochastic. 

It is important to distinguish between two different 
kinds of stochasticity. Wave-particle interaction in pres- 
ence of of uncorrelated small amplitude electromagnetic 
waves or plasma turbulence is called extrinsically diffu- 
sive [ll|- In this case the regular phase space structure 
for a charged particle interacting resonantly with an elec- 
tromagnetic wave is perturbed by neighboring uncorre- 
lated waves. This leads to extrinsic stochasticity and 
diffusive behavior. On the other hand nonlinear systems, 
such as particle interacting resonantly with a large am- 
plitude obliquely propagating (with respect to Bq) elec- 
tromagnetic plasma wave, can exhibit intrinsic stochas- 
ticity. Indeed, when the wave amplitude is sufficiently 
large, the resonances at the gyrofrequency harmonics are 
sufficiently broadened that they overlap with adjacent 
primary resonances. Therefore particles interacting even 
with a single monochromatic wave may exhibit intrinsi- 
cally stochastic and diffusive behavior [l^. This is the 
regime of nonlinear diffusion and irreversible chaotic mix- 
ing of orbits. 

Because one of the main hypothesis of quasilinear the- 
ory is that particles dynamics is adequately modeled by 
their unperturbed trajectories, the quasilinear timescale 
Tc must be much smaller than the timescale for the onset 
of nonlinear orbit effects Tni [cf. [l^ : 

Tc < T„; — , (5) 

where uji, is the bounce frequency. This means that the 
turbulent spectrum should be broad enough so that the 
typical timescale for a charged particle to interact with 
a resonant wave-packet would be much less than its typ- 
ical bounce time, — 27r/a;fc, in a monochromatic wave 
at the characteristic wavenumber and frequency of the 
wave-packet. The bounce time, Tf,, for a particle in reso- 
nance with an electromagnetic wave is proportional to its 
oscillation period in the pseudo-potential well governing 
the resonant wave-particle interaction [121] . This interac- 
tion can be approximated by a Hamiltonian pendulum in 
the vicinity of the resonance point. 

Particles in resonance with a single finite amplitude 
fluctuation undergo a finite amplitude nonlinear oscilla- 
tion. This is the so-called trapping width, At;|| , given 
by the half peak-to-peak difference in the particle ve- 
locity parallel component. The trapping width and the 
bounce frequency for a nonrelativistic particle interact- 
ing resonantly with an electror nag netic wave are given 
by Equations (5a)-(5c) of Ref. [l^. These approximate 
expressions for Aw|| and Wfc yield considerable physical 
insight into the diffusion process [l^ when used in con- 
junction with the quasilinear diffusion coefficient. 
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Arbitrary length scale 
Alfven speed 
Unit transit time 
Magnetic field 
Electric field 



A 

VA 

TA = X/VA 

Bo = V47rp va 

En = {va/c)Bo = V\y/Wp/C 



TABLE I. Characteristic physical quantities. 



III. MAGNETIC MOMENT TRAPPING WIDTH 

From the trapping width, Av\\ , and bounce frequency, 
Wf), computed by Ref. 16] for the case of a circularly po- 
larized electromagnetic wave (see Appendix), it is possi- 
ble to derive the pitch angle trapping half width as: 



Aa 



Avn 



Bo 



1/2 



(6) 



As magnetic moment /x is related to a by Eq. (|4]), we can 
write the trapping width for the magnetic moment as: 



SB 

Bn 



-,1/2 



(7) 



These expressions apply to a circularly polarized wave. 
From Eq. ^ we expect that /i continues to be a good 
adiabatic invariant when resonances are not present or 
when particle interacts with extremely small amplitude 
waves. 



IV. MODEL AND GOVERNING EQUATIONS 

We investigate magnetic moment behavior first during 
the resonant interaction between one ion and a circularly 
polarized magnetic wave, then when resonance overlap- 
ping occurs and finally during the interaction between 
a distribution of particles and a broad-band turbulent 
spectrum. Because some of our normalization quantities 
are expressed in terms of typical time and length scales of 
the turbulence slab model p, [l3| , we first give a general 
summary of the slab model. 

For the general one dimensional (ID) slab descrip- 
tion, turbulence is made up of a sum of right and left 
handed circularly polarized nondispersive plane Alfven 
waves propagating in the parallel direction. The mag- 
netic field fluctuations are perpendicular to both the wave 
vector and the mean field. The fields are assumed to be 
magnetostatic. This amounts to the auxiliary assump- 
tion that the average particle speed is well in excess of the 
phase speed of the underlying linear wave mode. We ig- 
nore nonlinear wave- wave couplings in the spirit of quasi- 
linear theory [see e.g.. [T8l - [20| . 

Considering Alfven waves propagating with w/fc = 
a;/fc|| ~ ivA, the magnetostatic approximation implies 
|v| ^ VA (strictly \v\\\ ^ va)- Since particle energy is 
conserved in a frame moving at the parallel component 



of the phase velocity of the wave (a;/fc|[), quasilinear the- 
ory [131 implies: 

(w|| — w/fc||)^ + = const. 

Because of the magnetostatic assumption, particle en- 
ergy is conserved, i.e., energy diffusion in forbidden and 
in velocity space the resonant interaction diffuses pitch 
angle and gyrophase only. Finally, we ignore all inter- 
particle correlations resulting from their mutual interac- 
tion through their microfields (e.g., Coulomb collisions, 
Debye shielding, and polarization). Furthermore the 
feedback of the particles on the macroscopic fields is ig- 
nored, i.e., we consider only test particles in prescribed 
macroscopic magnetostatic fields. By virtue of the in- 
equality U|| 3> VA, the turbulent electric field of the or- 
der {6B/Bo)vaBo is negligible compared to the motional 
electric field of the particle, v\\Bo. 

The dispersionless hypothesis rules out phase mix- 
ing and, hence, phase decorrelation due to this pro- 
cess. Consequently the only way for a particle to see 
a "wavepacket" phase-decorrelate is to traverse an auto- 
correlation length of the turbulence [2l[ . The autocorre- 
lation time in this case is given by 



1 



1 



|A(w-A:||W||)| |w||A/c||| {v^ 



vul 



(8) 



where Ac is the turbulence correlation length. 

The behavior of a test particle is described by its time 
dependent position r(t) and three-dimensional velocity 
v(t), that are advanced according to dv/dt — v and the 
Lorentz force equation: 



E- 



(9) 



In order to render the equations non-dimensional, we 
use the characteristic quantities listed in Table HI where 
TA is the Alfven crossing time, va is the Alfven veloc- 
ity, A = is the turbulence coherence length related 
to the turbulence correlation length Ac (Ac = 0.7471^ for 
our particular slab configuration ) . For the static case 
also the light speed may be used as a characteristic quan- 
tity . The introduction of an Alfven speed in our test 
particle model, where the waves are treated as static, may 
appear rather artificial. However, the magnetostatic as- 
sumption is valid here provided that jwn | ^ va and we 
introduce va in anticipation of future work where we will 
drop the magnetostatic hypothesis. 

With our choice for the characteristic quantities (Ta- 
ble IJ) the dimensionless equations of motion of our 
charged test particles are given by: 



rfr 
(ft 



(E-f V X B) 



(10) 
(11) 



Here /3 = Qta [cf. a parameter in Ref. l23| couples parti- 
cle and field spatial and temporal scales and provides a 
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particularly useful means to relate our numerical exper- 
iments to space and astrophysical plasmas. In general 
in a turbulent collisionless plasma the bandwidth of the 
inertial range fluctuations may extend from large fluctua- 
tions at the correlation scale, Ac, to small fluctuations at 
the ion inertial scale. In this case /? ^ I and the turbu- 
lent time-scales are much slower than the typical particle 
gyroradius [25]. 

The resonant condition for the static case in terms of 
P is given by 



g b = 0.001 



8 b = 0.01 



gb = 0.1 



n(3 



n(3 



(12) 



Time is advanced through a fourth-order Runge-Kutta 
integration method with an adaptive time-step [pp. 708- 
716 of Ref.'26'l. 



V. NUMERICAL SIMULATIONS 

Particles are loaded randomly in space at < = 
throughout a one-dimensional simulation box of length 
L. The fields are described in the following sections. In 
spherical coordinates, with the polar axis along the z- 
direction parallel to the mean magnetic field of strength 
Bq, particle velocity components are: 



V sm a cos < 



t;sin0sin0 Vz = vcos6 (13) 



Particles initial velocities are randomly distributed in the 
gyrophase (j) between [0 : 27r], while the velocity magni- 
tude V and pitch angle 6 are determined by the particular 
numerical experiment. 

Typical particle velocities used in our simulations are 
IQvA and IQQva, satisfying the magnetostatic constraint. 
In our analysis magnetic moments are expressed in units 
of the characteristic quantity = v'^ /Bq. We also define 
Sb = SB/Bo. 

The statistic analysis of particle magnetic moment in- 
volves averaging trajectories over the particle gyroperiod 
Tg = 27r/f2. For each simulation we compute the effective 
number of gyroperiods N^- that particles complete in a 
given magnetic field configuration as: 



* dt eB{t) 
27r mc 



(14) 



where B{t) is the intensity of the total magnetic field. 
When 5h ^ 1, B{t) ~ Bq] however increasing 5h toward 
unity the waves contribution to the strength of the total 
magnetic field B{t) is not negligible. 



V.l. Single wave 

We start studying the ion motion in presence of a con- 
stant magnetic field Bg and a perpendicular left-handed 
circularly polarized wave with 

B — 5Bx cos{k()z) — 6By sin (fco^:) e^^ -I- Bq e^, (15) 




FIG. 1. Gyroresonant interaction between a circularly po- 
larized wave and a particle with i; = IOOha and a = 1/8: 
cosine of pitch angle a (top row), particle magnetic moment 
fj, (middle row) and parallel component of the induced electric 
field bottom (row). Different columns correspond to different 
wave amplitude: Sb — 0.001 (first column), Sb = 0.01 (sec- 
ond column), Sb = 0.1 (third column) and 5b — 1.0 (fourth 
column) . 



where SBj. and SBy are the amplitudes of the wave and 
fco is the wavevector. We assume SBx = SBy = SB for 
the rms average values. In these simulations /3 = 10"^, 
V ^ lOOvA and a = 0.125 {0 = 82°). 

We follow the test-particles until they complete N^^ — 
100 gyroperiods. For the resonance condition, Eq. (fRt . 
we set ko = 80/A. Particles injected with a pitch angle 
cosine different to a = 0.125 will not be in resonance 
with this wave, exhibiting a different behavior. For a 
direct comparison we also inject non-resonant particles, 
i.e., with a = 0.5 {9 = 60°). 

Figure [T] shows the time evolution of the cosine of pitch 
angle a, particle magnetic moment /z, and the parallel 
component of the induced electric field , for a resonant 
particle (a = 0.125). Different columns corresponds to 
different values of the wave amplitude: Sb ~ 0.001 (first 
column), Sb = 0.01 (second column), Sb = 0.1 (third 
column) and Sb — 1.0 (fourth column). 

When the parallel component of the induced electric 
field is almost constant and equal to ^ -'V±Sb, the 
resonant interaction produces variations that are secular 
over a gyroperiod. However an oscillation occurs over a 
longer time, the bounce period T}, — 2tt/uji, (where is 
the bounce frequency discussed in Section [TT|. This is the 
typical timescale over which the velocity, and hence the 
particle trajectory, exhibits significant deviations from 
the linear = const and v± = const case. 

In Section IIIII we derived the analytical expression 
for the half trapping-width of magnetic moment for a 
particle interacting with a left or right handed circu- 
larly polarized wave (see Eq. [7]). We now compute the 
values of the half peak-to-peak difference in a and /z, 
Aa = {amax - Q;mm)/2 and A/u, = {umax - IJ'min)/'2, for 



5 



Sb 










0.001 


0.022 


0.02 


0.0056 


0.0055 


0.01 


0.07 


0.075 


0.0176 


0.02 


0.1 


0.2227 


0.2 


0.0556 


0.055 


1.0 


0.704 


0.6 


0.176 


0.175 



TABLE II. Trapping width values for a and /i: comparison 
between theoretical (subscript th) and numerical (subscript 
sim) values. 



the resonant interaction simulations. These values and 
those obtained from the theoretical expressions (|6l)-(l7]) 
are listed in Table and are in good agreement, con- 
firming the validity of equations (|ni)-(IZl) and reinforcing 
the intuitively idea that magnetic moment and pitch an- 
gle behaviors are strictly related. 

To compare resonant and non-resonant dynamics, we 
show in Figure [5] the time evolution of cosine of pitch 
angle a (first row), magnetic moment /i (third row), and 
their distribution functions f{a) (second row) and /(/i) 
(fourth row) at the end of the simulation, for a reso- 
nant particle with a = 0.125 (left column), and a non- 
resonant one with a = 0.5 (right column). In contrast 
with the resonant case in which a and fi exhibit well- 
known secular variations with typical period equal to t^, 
the a and /i profiles for a non-resonant particle show 
a regular oscillating behavior, a distinctive signature of 
regular particle motion. The values of the half peak-to- 
peak difference in a and /i obtained from the simulation 
are Aofsim = 0.0025 and Afigim = 0.003. These are 
smaller than the theoretical values computed from equa- 
tions (in|)-(I7]) with 6b = 0.01 and a = 0.5, for which we 
obtain Aath = 0.1316 and Afith = 0.1316. 

The distribution functions /(a) and /(/x) (Figure ^ 
for a resonant particle are more spread in a and /i and 
are centered around their initial values a = 0.125 and 
/i = 0.98. In the non-resonant case, /(/i) remains peaked 
at its initial value, i.e., its magnetic moment is constant 
during particle motion. The spread in a of its distri- 
bution is ^ 10%, small compared to the resonant case 
spreading of ^ 40%. 

Figure [3] shows the distribution functions, /(a) and 
f{fi), at the end of the simulation for 1000 resonant and 
non-resonant particles injected in the simulation box with 
random positions and phases. For non-resonant particles 
(right column) the distributions remain peaked around 
their initial values a = 0.5 and /i//in = 0.75 with very 
little spreading. For the resonant particles (left column) 
/(a) acquires a Gaussian shape centered around its ini- 
tial value a = 0.125. Furthermore it spreads of ~ 0.1, 
comparable to the trapping width for the single particle 
2Aa = 0.014 (Figure [5]). The magnetic moment distri- 
bution for the resonant case has a characteristic shape 
found for /i in the parameter range in which pitch angle 
exhibits a Gaussian distribution and the density distri- 
bution function is still isotropic (particle free-streaming 
regime). As for the pitch angle, the spread in the mag- 
netic moment distribution of 0.03 is comparable to the 
trapping width for the single particle 2A/x = 0.00352 (see 
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FIG. 2. Time evolution of cosine of pitch angle a (first row) 
and its distribution function /(a) (second row), time evo- 
lution of magnetic moment fj, (third row) and its distribu- 
tion function /(/i) (fourth row) of resonant {a = 0.125, left 
column) and non-resonant particle (a = 0.5, right column). 

V = 100 VA. 



Eq.[7]and Figured]). 



V.2. Overlapping resonances 

In order to understand the effect of overlapping reso- 
nances on particle magnetic moment, we perform a nu- 
merical experiment with four different particles in the 
simulation box with random initial positions, same ini- 
tial velocity v = 100 va, but different values for pitch 
angle cosine: ai = 1/2, a2 — 1/4, as = 1/8, 04 — 1/32. 
For (3 = 10^, making use of the resonance condition for 
the static case [Eq. ([T^ ]. the cyclotron resonances n = 1 
for the different values of a are expected for fciA = 20, 
fczA = 40, fcsA = 80, and = 320. 

The total magnetic field is given by: 



B = BqBz -\- Sbcos[kiZ -\- (f>i]ex — Sbsm[kiZ -t- 4>i]^y, 

i=l i=l 

(16) 

where the 4>i are random phases. Taking into account 
resonance broadening effects, all particles with parallel 
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FIG. 3. f(a) (first row) and f{n) (second row) at the end of 
tlie simulation for an initial distribution of 1000 resonant (left 
column) and non-resonant (right column) particles randomly 
distributed in the simulation box. Sb = 0.01. 



Sb 


Q = 1/2 


Q = 1/4 


a = 1/8 


Q = 1/32 


0.001 


4.16 


3.1 


2.227 


1.3 


0.01 


13.1 


9.85 


7.042 


3.583 


0.1 


41 


31.1 


22.27 


11.33 


1.0 


131 


98.3 


70.42 


35.83 



TABLE III. Values of Auy for a =1/2, 1/4, 1/8, and 1/32 
resonances at different Sb. 
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8b=1.0 
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FIG. 4. Transition from non-overlapping to overlapping reso- 
nances: a (left column) and /j, (right column) profiles varying 
the waves amplitude: Sb = 0.001 (first row), Sb — 0.01 (second 
row), Sb = 0.1 (third row), Sb — 1.0 (fourth row). 



velocities in the range 

D|l — Aw|| < I'll < -I- Awy (17) 

can potentially resonate with a wave, whose wave number 
is fc|| = il/v\\. As found by Ref. [g^l, the direct evidence 
of resonances overlapping is the disappearance of con- 
stants of motion, i.e., the onset of stochasticity in the 
Hamiltonian formalism. We make simulations with four 
different waves amplitudes Sb = 0.001, 0.01, 0.1, and 1.0. 
The values of the trapping half- widths Ai;|| computed 
for the different pitch angles with Eq. (jA.SP are listed in 
Table Hill for the different 6b considered. 

Figure |4] shows time histories of pitch angle cosine a 
(left column) and magnetic moment fi (right column) 
profiles for various 6b. Again similar behavior is seen for 
a and fi. For the smallest wave amplitude, 6b = 0.001, it 
is possible to recognize very well the four different reso- 
nances in the profiles of a and /i. For 6b — 0.01 the res- 
onance at aa = 1 /8 is overlapping with the resonance at 
a4 = 1/32. Indeed, the initial parallel velocity of the par- 
ticle injected at the smallest pitch angle, W|| 4 = 3.125uyi, 
lies in the range of velocities [see Eq. (fT7| ] in possible 



resonance with fc|| — k^. For higher wave amplitudes, 
6b — 0.1 and 6b = 1.0, the condition (fTT]) is satisfied 
by all particles velocities. Stochasticity arises and the 
different resonances are indistinguishable. 

The distribution functions /(a), fifJ-) and f{6z) (where 
6z — z ~ zq is the displacement along z relative to the 
particle initial position zq) after lOOrg (Figure [5]) exhibit 
similar characteristics. 

For 6b = 0.001, /(a) and /(^) are peaked in corre- 
spondence of their four initial values because of the good 
resonances separation. f{6z) shows that the particles 
are simply free-streaming in the parallel direction and, 
depending on their initial parallel velocity, they cover 
shorter or longer distances along z. 

For 6b = 0.01, f{a) spreads around its initial four 
peaks because particle interact resonantly with waves of 
larger amplitude, and resonances overlap for a < 1/4, as 
discussed previously. Similar effects are shown also by 
/(/i), confirming that for small 6b the resonant interac- 
tion affects magnetic moment and pitch angle in similar 
ways. 

While for 6b = 0.01 particles continue to free-stream 
in the z-direction, different profiles for f(5z) appear for 
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FIG. 5. Transition from non-overlapping to overlapping res- 
onances with varying wave amplitude: Sb = 0.001 (first row), 
5b = 0.01 (second row), Sb = 0.1 (third row), Sb = 1.0 (fourth 
row). The distribution functions f{a) (left column), /(/i) 
(central column) and f{Sz) (right column) are averaged over 
time. 



Sb = 0.1. Pitch angle distribution begins to isotropize 
and magnetic moment exhibits a one-sided long tail dis- 
tribution extending toward smaller /x. This behavior is 
similar to the regime found previously in the single wave 
experiment when /(a) is nearly isotropic, f{Sz) still indi- 
cates particles free-streaming, and the magnetic moment 
distribution displays a long tail. 

For Sb = 1.0, by < = IOOtc, the pitch angle cosine distri- 
bution f{a) has become completely isotropic, while f{Sz) 
approaches a gaussian distribution indicative of spatial 
diffusion. In this regime /(/i) loses its long-tail and starts 
to acquire a gaussian shape. In that way we have identi- 
fied three distinct regimes of statistical magnetic moment 
behavior with increasing degree of turbulence. 



VI. SLAB SPECTRUM 

In this section we present the results of our numerical 
simulations of test-particles in presence of a broad-band 
slab spectrum [see Eq. and Figure [S] . We have per- 
formed simulations for different particles velocities and 
amplitude of the magnetic field fluctuations. 

Simulations use a unidimensional computational box 
of length L = 10000 Iz {Iz = 1 is the coherence scale for 
the slab spectrum) with Nz = = 268,435,456 grid 
points. The magnetic field in physical space is generated 
from a spectrum P{k) in Fourier space, via inverse fast 
Fourier transform (FFT). The turbulent magnetic field 
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FIG. 6. Power spectrum of the turbulent magnetic, k is 
normalized to the coherence length h. 



is given by: 



B{z) = Baez+6B{z), 



(18) 



with (5B(z) = 6Bx{z) ex + SBy{z) By and the solenoidality 
condition is identically satisfied. 

The modes of the magnetic field components in k-space 
are given by: 

5B,{K) = [P(fc„)]i/2e^*" 
5By{K) = [P(fc„)]i/2e^*" 

where fc„ = 2Tm/L and <&„ and are random phases. 
The slab spectrum P{k) is given by: 

{Cslab[l + (fcn^^)^]"^/^ for K < kdiss 
[^dissik^) ' for fc„ > fcdiss 

where Csiab ~ 2XcSb1 ^^^^ is a constant specific to this 
form of the slab model, Sb^ ^.^^^ is the mean square 
fiuctuation, kdiss is the dissipation range wavenumber, 
Cdiss = Csiab[^ + (kdisslz)'^]'^^^ is the constant for the 
dissipation range (set by the continuity of the spectrum 
P{k) at kdiss)- The vectors of Fourier coefficients are 
zero-padded for N„iax + ^ l£ n < Nz providing an 
extra level of smoothness to the fields by an effective 
trigonometric interpolation. In all the simulations we 
use Nraax = 6.7 X 10^ and a simple linear interpolation 
to compute the fields at the test particle position. 

The resulting spectrum is shown in Figure [6] Several 
important scales are present in the system. They are 
labeled as kmin, h^, kdiss, kmax and kN,- The discrete 
wavenumbers are obtained through fc„ = 21m /L as: 



Nk = ^k^ 1600 A:. 
2tt 



(20) 



We summarize the values for k and iV^ used in our 
simulations in Table IWl where: 

- kmin = 2'k/L is the minimum wave vector of the 
spectrum, corresponding to Nk = Nkmin = 1- 
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TABLE IV. Characteristic scales in the spectrum. 

Wavenumber index Wavenumber value 

A^fe™» = 1 fc„,„ = 6.28 X 10-* 
Nki, = 10^ k,^ = 27V 

A^fe^^^ = 1.6 X lO'' kaiss = 10^ 

^feMAx = 6.7 X 10^ fcMAx = 4.2 X 10* 

^fcMAxi = 1-3 X 10^ kMAxi = 8.4 X 10-* 

- ki^ = 2n/lz = 27r is the wave vector that marks 
the beginning of the inertial range. Three decades 
of energy containing range from kmin to fc;^ ensure 
turbulence homogeneity. 1^ or Ac = 0.747^2 [see 
Ref. [13 correspond to the typical lengths scales 
over which the particles attain diffusive behavior 
of the pitch angle. Three decades of inertial range 
with P{k) cx k~^/^ well represent solar wind condi- 
tions. 

- kdiss is the wave vector corresponding to the be- 
ginning of the dissipation range. In our model, the 
spectrum extends beyond kdiss with P(k) oc k~'^/^ . 

- A t two d ecades higher wavenumber, kMAX = 
^ffnjjmlkdiss determines the end of the dissipation 
range. 

- Extending for two decades beyond kMAX = 4.2 x 
lO**, the spectrum includes zero-padding up to 
kMAxi = 8.4 X 10^. 

- Another important scale, not labeled in Figure [S] 
because it depends on test-particle velocity, is the 
wave vector corresponding to z^ax — vTtot, the 
distance covered by a charged test particle mov- 
ing at speed v in the simulation running time Ttot- 
To avoid periodicity effects it is important that the 
box length L is large enough so that particles tra- 
jectories are limited to a small fraction of the full 
length, i.e., L > Zmax or kmin < ^/zmax- Period- 
icity might indeed give rise to artificial field lines 
diffusion. 

We fix the value of the /? parameter equal to 10^. This 
corresponds approximately to observed solar wind turbu- 
lence properties at 1 AU, as follows: 

\mJ c c An' 

where Ac is the turbulence correlation length ujpi = 
(47rnoi(7i^/mi)^/^ is the ion plasma frequency (qi and 
rui are respectively the ion charge and mass) and Xn = 

1/2 

c/ujpi = (yC^mi/ AimiC^^ is the ion inertial length. For 
Ui — He then \ii — {mi/me)^^'^ Pie- Because the solar 
wind density at 1 AU is approximately n ^ (Ij 10) cm""^ 
on average A^^ ~ 1000 km. At the same distance the 
turbulence correlation length Ac is approximately 10^ km 
[H] and /? ~ 10^. 

Typically 1000 particles are injected in the simulation 
with initial random positions. Particles are loaded from 



TABLE V. Typical values used in the simulations. 



V\va\ 


rL[h] 


hires 


e 


tc[TA\ 


10 
100 


10-2 


8 X 10^ 
8 X 10^ 


1.33 X 10-^ 
1.33 X 10-2 


0.0747 
0.00747 



a cold ring beam [see equations (|13p ] distribution with 
constant velocity magnitude, sin0 is set equal to (1 — 
Q!q)^/^, where uq is the initial pitch angle cosine respect 
to the background field Bq. The initial gyrophase is 
chosen randomly. For all the simulations ao = 0.125 
(6*- 82°). 

From the previous section, we know that the behavior 
of magnetic moment is correlated to pitch angle behavior 
for a low level of magnetic fluctuation {5b = 0.001, 0.01). 
Pitch angle and magnetic moment exhibit Gaussian dis- 
tribution functions typical of normal diffusion processes. 
Increasing the turbulence level, pitch angle distribution 
approaches isotropization and a transient regime is ob- 
served with the magnetic moment starting to be influ- 
enced by the onset of spatial parallel diffusion. When 
f{a) completely isotropizes, spatial diffusion sets in and 
/(/i) behavior is closely related to the sampling of the 
varying magnetic field strength associated with that spa- 
tial diffusion. 

From quasilinear theory we know that velocity and 
real space diffusion occur at two different time scales. 
Typically, velocity space diffusion takes place with the 
time scale Tc = \c/v shorter than the typical time scale 
at which parallel diffusion occurs T|| = A||/w, where 
A|| — 3-D||/t; is the parallel mean free path. For this 
reason we follow test particles in the simulation box for 
a time T > Tc, typically with T = 20tc. Particles param- 
eters used in the simulations are listed in Table |Vl 

An important parameter in the description of energetic 
test particles is e = r^/Xc, which is sometimes called the 
dimensionless particle rigidity. It can be related to the 
bend-over wavenumber of the turbulence, kbo = 1/Ac, 
and the minimum resonant wavenumber, fcj„i„ = I/^'l, 
as e = kbolKnin- For example when ^ Ac particles 
experience all possible fc-modes in few gyroperiods res- 
onating with the energy containing scale (fc^j„ ^ fcfjo)- 
For lower energies the test particles resonate in the iner- 
tial range. Those with v = 1Qva will resonate at the end 
of the inertial range (1/ri in Fig. 15]), while those with 
V — \QQva at the middle of the inertial range (l/r2 in 
Fig.ini). Furthermore, as explained previously, the condi- 
tion kmin ^ Zmax is necessary to avoid artificial effects 
in particle transport associated with periodicity of the 
magnetic field. 

Figure [7] shows /(a) (left column), /(/i) (central col- 
umn) and f{Sz) (right column) for a distribution of parti- 
cles moving with an initial velocity v = IOva in presence 
of the slab spectrum [Eq. (|19p , Figure|B] . All the distribu- 
tion functions are computed at the end of the simulation, 
i.e., after 20 Tc. The blue line and the green line indicate 
the initial value and the mean value of each distribution. 
As particles are injected at different positions, it is con- 
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FIG. 7. Distribution functions of cosine of pitch angle f{a) (left column), magnetic moment /(/i) (central column) and particle 
displacements relative to the initial position f{5z) (right column) at 20 Tc for different waves amplitude: 5b = 0.001 (first row), 
Sb = 0.01 (second row), 5b — 0.1 (third row) and 5b = 1.0 (fourth row). Particle parameters at injection: v — Wva and 
ao = 0.125. 
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FIG. 8. Statistics for v = Wva- Variances for cosine of pitch angle a (left plot), magnetic moment /i (center plot) and particle 
displacements relative to the initial position 5z (right plot) at different values of 5b: 5b = 0.001 (black line), 5b = 0.01 (purple 
line), 5b = 0.1 (red line) and 5b — 1.0 (blue line). The s values indicate the different slopes for the scalings ((Aa)^), {(A/i)^) 
and {{Azf) oc t" . The variances are fitted with the gray lines: dotted line is s = 2, dotted-dashed line s = 1, dashed line 
s — 0.8, three dotted-dashed line s = 0.7. 
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venient to define the quantity 8z = z[j) — z(0) [j is a 
temporal index). In this way it is possible to take out 
from the distribution function /((5z) both the drift effect 
(vut) and particle diffusion relative to their own posi- 
tions (Azj). The general expression for the z position of 
the z-th particle is given by 

z^ = Zi{0) + VDt + Az, = Zi{0) + Szi. (22) 

The primary diagnostic for studying particle diffusion 
is the variance {t) cx t'' of particles cosine of pitch an- 
gle, magnetic moment and position parallel to the mean 
field direction. Figure [8] illustrates the time evolution of 
the variances, ((Aa)^) (left figure), ((A/i) ) (central fig- 
ure) and ((Az)^) (left figure) for a particles distribution 
moving with initial velocity equal to 10 v a- Different col- 
ors correspond to different 6b values: 6b = 0.001 black 
line, 6b = 0.01 purple line, 6b = 0.1 red hne and 6b = 1.0 
blue line. The variances are fitted with the gray lines: 
the dotted line is used for s = 2, the dotted-dashed line 
for s = 1, the dashed line for s — 0.8 and the three 
dotted-dashed line for s = 0.7. 

For 6b = 0.001, a and /i display Gaussian distributions 
while particles free-stream in the z-direction. Particles 
that cover greater distance in z are more scattered in 
pitch angle and consequently in fi. Figure [8] shows su- 
perdiffusive behavior (black line, s — 2) with particles 
free streaming along z, and later^ variance characteris- 
tic of normal diffusion with ((Aa) ) and ((A/i) ) scaling 
(X t. 

For 6b = 0.01, particles cover only one side of the a 
hemisphere continuing to travel along z (purple line in 
Figure [8|). This is the transient regime already observed 
in Figure [5] when /(/i) exhibits a one-sided long tail dis- 
tribution toward smaller /i. 

For 6b = 0.1, pitch angle distribution becomes com- 
pletely isotropic and spatial diffusion sets in, as shown 
by the slope s = 1 of ((Az)^) in Figure [5] (red line) at the 
end of the simulation. The deviation from purely free- 
streaming or ballistic behavior means that, while the sys- 
tem has not become fully diffusive along the mean field 
direction, there are signs that diffusive processes in ve- 
locity space are beginning to diminish the free-streaming. 
Although /(yu) still exhibits a long-tail, the influence of 
spatial diffusion starts to appear. The well-pronounced 
peak observed in f{fi) for 6b = 0.01 is substantially re- 
duced and the mean value of magnetic moment decreases. 
Moreover fi displays subdiffusive behavior up to 0.02tc. 
After this time particles diffuse in space and ((A/i)^) at- 
tains a plateau. The Gaussian shape is not reached yet 
probably because spatial diffusion is just at the begin- 
ning. 

For 6b = 0.5 [see Figure ®] and 6b = 1.0, f{a) is 
isotropic, particle motion is conyletely diffusive in real 
space [as the slope s = 1 in ((Az) ) in Figure[8](blue hne) 
shows], and f{ii) behavior is closely related to the sam- 
pling of varying magnetic field strength associated with 
that spatial diffusion, displaying a Gaussian distribution 
centered at the middle of /i-space. 

From a more detailed analysis of the case 6b = 0.5 



[Figure ([9])] we notice that magnetic moment variance 
(first figure) scales according to ((A/i) ) cx t°-^'^ (red- 
dashed line) up to 0.002tc and after 0.005tc a plateau 
is attained (blue-dashed line). Instead particles motion 
(see time evolution of ((Az)^), second plot) becomes fully 
diffusive (blue-dashed line) only after O.OOTtc. Magnetic 
moment distribution /(/i) in the first part of the evolu- 
tion (third plot) is in the transient regime characterized 
by the long tail. In contrast, when particles diffuse in real 
space (fourth plot), f{fJ.) reacquires the gaussian profile. 
Thus the final stage of magnetic moment variance evolu- 
tion, i.e. the formation of the plateau, can be considered 
as a precursor for the onset of the parallel diffusion of par- 
ticles in space. Of course this effect is present in pitch 
angle variance too, but in addition in /i behavior we have 
a direct signature of the onset of the spatial diffusion, 
that is the reappearance of the gaussian shape in the dis- 
tribution function, while pitch angle distribution remains 
completely isotropic. 

Thus these transitions in magnetic moment behavior 
are related not just to the variation of the turbulence 
level, but also to the different time scale at which mag- 
netic moment conservation is studied. 

The magnetic moment distribution functions /(/i) and 
variances ((A/i) ) in the case v = 100 (not shown) ex- 
hibit the same features observed for v — IOva- However, 
increasing particle speed the total number of gyroperiods, 
Nrg, performed by each particle decreases; as a conse- 
quence, faster particles sample less variation in magnetic 
field strength. This leads to a slower spatial diffusion, 
i.e., for lOOvA spatial diffusion occurs on a time scale 
longer than 20rc. 

For V = lOOvA we show in Figure (|10p magnetic mo- 
ment standard deviation a^/^min (blue triangle) and the 
changes in its mean A/i//ii„ = {fl — ^Xm)/ p^in versus 6b 
after 20rc. As 6b increases toward unity the changes in 
magnetic moment distribution start to increase faster. 



VII. CONCLUSIONS 

In this paper we have investigated the conservation of 
charged particle magnetic moment in presence of turbu- 
lent magnetic fields. For slow spatial and temporal varia- 
tions of the magnetic field respect to the particle gyrora- 
dius and gyroperiod, the magnetic moment /i is an adia- 
batic invariant of the particle motion. Non-conservation 
of magnetic moment can influence particle acceleration 
and have considerable implications in many astrophysical 
problems such as coronal heating, cosmic rays transport 
and temperature anisotropics in the solar wind. These 
applications motivate the present basic study of the de- 
gree to which magnetic moments are conserved in increas- 
ingly complex models of one dimensional spectra. While 
all the models considered here have been very oversim- 
plified relative to the spectra observed for example in the 
solar wind [13, [13 O'" simulations of MHD turbulence 
[30I Isij . the present study is intended to contribute to 
the basic understanding of the conditions for the onset 
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of magnetic moment nonconservation. We point the in- 
terested reader also to a recent study by Ref. [s^ that 
addresses this issue from a somewhat different perspec- 
tive. 

In order to reproduce and extend some of the result 
15|, we started to study the resonant 



obtained by Ref. 
interaction between ions and a single parallel propagat- 
ing electromagnetic wave (see Section |V.1[) . Using the 
specialized expression for the trapping width Aw|| found 
by Ref. [l^ in the case of a single circularly polarized 
wave, we have been able to write a similar expression for 
magnetic moment (see Eq. [T]). In presence of a single 
finite amplitude fluctuation the magnetic moment of a 
resonant particle undergoes a finite amplitude nonlinear 
oscillation too. We have performed several simulations 
changing both particle velocity and the amplitude of the 
wave. For each of them we compare the values of A/i and 
Au|| with those obtained using our specialized expression 
and they are in good agreement. 

We designed a particular experiment to study the ef- 

From 



fects of resonances overlapping (see Section IV.2|) 
the analysis of the distribution functions of particles pitch 
angle, f{cy), magnetic moment, /(/i), and z-position, 
/(z), we distinguish three different regimes. First, 
for a low level of magnetic fluctuation, i.e., SB/Bq = 
0.001, 0.01, the magnetic moment distribution half-width 
is directly related to pitch angle distribution. Second 



for SB/Bq = 0.1 stochasticity arises as a consequence 
of overlapping resonances and its effect on pitch angle 
is isotropization of the distribution function. This is a 
transient regime during which magnetic moment exhibits 
a one-sided long-tail distribution and starts to be influ- 
enced by the onset of spatial parallel diffusion. Finally, 
when /(a) completely isotropizes spatial diffusion sets in 
( SB/Bq = 0.1), /(^) behavior is closely related to the 
sampling of varying magnetic field strength associated 
with that spatial diffusion. 

Other studies regarding particles interaction with two 
electromagnetic waves as well as a flat turbulent spec- 
trum (not shown) were also conducted and they con- 
firmed this general picture. 

Motivated by these results we studied the behavior of 
many particles interacting with a broad-band slab spec- 
trum, generated in order to mimic some of the major fea- 
tures of the solar wind (see Figure ^ : (a) three decades 
of the energy containing scale ensure turbulence homo- 
geneity, (b) three decades of inertial range well-reproduce 
the observations and (c) two decades of dissipation range 
enable us to cross the "amin barrier" related with the 
"resonance gap" predicted by quasilinear theory [2l|, [s^] . 
After that there are almost other two decades of zero- 
padding, important for the trigonometric interpolations 
and for the smoothness of the field. This is implemented 
using a numerical grid with = 2^^ — 268,435,456 
points corresponding to 134, 217, 728 wavevectors for the 
spectrum. Apart from the obvious limitation that this 
spectrum is purely one dimensional, it is constructed to 
correspond roughly to features of solar wind spectra ob- 
served by single spacecraft, where the fully three dimen- 
sional spectrum is in effect reduced to a one dimensional 
form. Information is lost in the process [see, e.g.. [29j. 

In order to gain insight on magnetic moment con- 
servation we have performed simulations changing 
both particles velocity, v = (10, 100) va, and the 
amplitude of magnetic field fluctuations SB/Bq = 
(0.001,0.01,0.1,0.5,1.0). Particles injected at different 
velocities start to resonate at different points of the spec- 
trum. We analyzed the distribution function [see Fig- 
ure [7] and the variance [see Figure (|5])] of pitch angle 
cosine a, magnetic moment /x and parallel position z. 
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From the experiment of resonances overlapping we 
know that the three different regimes of /i statistical be- 
havior are related with other two effects; diffusion in 
velocity space and spatial parallel diffusion. These take 
place at different characteristic times, Tc and T|| respec- 
tively. In order to investigate the effects of both pro- 
cesses on magnetic moment distributions, we followed 
test-particles in the simulation box for times T > t,.. 

For a low level of magnetic fluctuations particles free- 
stream in the z-direction while a and /i exhibit gaussian 
distributions around their initial values. For 6B/B0 — 
0.01 particles cover completely one side of the a hemi- 
sphere continuing to stream freely along z. This is the 
transient regime during which /(/i) exhibits a one-sided 
long tail distribution in the direction of smaller fi that 
appears to be a typical feature of magnetic moment dis- 
tribution. During this transient regime the distribution 
of particles nearly conserves its magnetic moment. In- 
creasing the value of SB/Bq spatial diffusion starts to 
take place, f{fJ.) recovers the typical Gaussian shape cen- 
tered in the middle of /i-space. These different regimes 
of magnetic moment statistical behavior are related not 
just to the variation of the turbulence level SB/Bq, but 
also to the different time scale at which magnetic moment 
conservation is studied [see Figure ©]. 

In spite of the limitations of the present approach the 
results presented here provide a basic view of how mag- 
netic moments are modified in simplified models, and in 
particular how magnetic moment changes are related to 
pitch angle changes and sampling of magnetic variations 
due to spatial diffusion. It is clear that additional study 
is required to understand more fully the influences of 
turbulence on magnetic moment statistics. For exam- 
ple realistic three dimensional models of the magnetic 
field turbulence, as well as incorporation of electric field 
fiuctuation effects, are expected to have significant ef- 
fects. It is also possible that nonGaussian features of 
magnetic field fluctuations, such as, are associated with 
intermittency effects, may also influence magnetic mo- 
ment changes, much as they influence spatial transport 
due to trapping and related influences [H, HI] ■ In this 
regard the present results, along with those of Ref. [30| . 
may be considered as baseline or minimal quantification 
of nonconservation of magnetic moments of a distribution 
of test particles in turbulence. Planned future studies will 
investigate quantitatively how additional realism in the 
modeling might produce even more significant departures 
from magnetic moment conservation. 
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Appendix: Derivation of trapping half width for a 
circularly polarized wave 

Using equations (5a) and (56) of [lH| it is possible to 
derive a simplified expression for the trapping half- width 
and the bounce frequency in the case of an Alfven static 
wave [? ]. For this particular case k± — and <j> = 0. 
We can rewrite equation (5 c) of [l^ as 

k 



2c 



£2 



£2 + -^^^1 ) Jn+l{k^p) 



;||K_L 
ck 



ei + £3 J7i{k±p) 



(A.l) 



with cos a = 1 and sin a = 0. Because k || Bq we 
can choose Bz = Bo/|Bo|, Gy is any arbitrary direc- 
tion perpendicular to and e^; = x e^. The vec- 
tor potential can be obtained from the magnetic field 
V X Bj^ BxBx + ByBy. In Fourier space V 
we have: 



ikz&z, so 



Ax - —By 



Ay - -Bx 



(A.2) 



Considering only a single circularly polarized wave in 
space, for the two different possible helicities we can 
write: 



B4 



(S±e±) exp [i{k\\: 



(A.3) 



where 

B± = ^{BxTiBy) and e± = -^(ex T ie^) (A.4) 

are respectively the complex amplitudes and the orthog- 
onal polarization unit vectors. The +(— ) polarization 
state is the positive (negative) helicity, i.e., the vector B 
is rotating counter-clockwise (clockwise). At first, let's 
consider only the left-handed polarized wave B+. As- 
suming B+ — v^f^iJe^"/^ we can write the x and y 
components of the wave magnetic field as 

Bx = SBexp [i(fc||Z - 7r/2)] 

By = SB exp {ik\i z) (A. 5) 

Inserting this two expressions into Eg. IA.2l we obtain: 



An, 



'II 

SB 



■W 



7r/2)] 



exp (ifc||z) 



Comparing the real parts of these equations with equa- 
tion (16) of [1^ we obtain an expression for the coeffi- 
cients Ai and A2 and for the normalized components of 
the wave polarization vector ei, £2 and £3: 



Ai =r/ 



SB 



A. 



SB 

k\\ 



where rj 



(A.' 
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TABLE VI. Wave polarization and resonance contribution to 
trapping width. 



Polarization 


V 


resonance n 


B+ left-handed 


1 parallel 


-1 


B- right-handed 


— 1 anti-parallel 


1 


B+ left-handed 


— 1 parallel 


1 


B- right-handed 


1 parallel 


-1 



_\q\rjSB \q\SB 



(A.7) 



Similarly, for a right-handed circularly polarized wave 
B_ we have: 



SB ^ SB , fc|| 
Ai = -rj—, A2 = — , where ?7 = 7^ 
fell K\\ fell 



ei 



\q\ri5B \q\5B 
£2 



mc^fc|| 



mc^kn 



£3 = 0. 



In case of a single circularly polarized wave propagating 
parallel (or antiparallel) to the magnetic field there is 
only one resonance present and particle motion is inte- 
grable [l^: indeed J„(0) ~ unless n = 0. Therefore 
depending on the polarization of the wave and on its di- 
rection of propagation only I = 1 01 I = —\ resonances 
contribute to the trapping width, as shown in Table IVTl 
Thus, considering equations (5a) and (56) of [1^1, Eg. lA.ll 
and Eg. rOlTOl with Jo(0) = 1, [? ] find a specialized 
formula for the trapping half width and bounce frequency 
applied to the case of a circularly polarized wave propa- 
gating parallel fc|| > and n = —1, or antiparallel, fc|| > 
and n = 1 to Bq: 



SB 



(l-a^/^H 

^0 

- (l-g^)!/^ SB^^^^ 
\a\ Bn 



1/2 



(A., 



if ^iit^ii > and zero otherwise, in which a = cos 9 is the 
cosine of pitch angle. Exactly the same set of equations 
holds for At'li'''^-* and uj^^'^^K However the condition for 
their being nonzero is reversed, i.e., fc||V|| > 0. We omit 
the superscripts (±1) because of this degeneracy. 
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